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A  method  is  used  in  which  the  discrete  vortex  model  is 
combined  with  boundary  layer  theory;  the  boundary  layer  separa¬ 
tion  and  the  rear  shear  layer  separation  are  simultaneously 
taken  into  account.  The  authors  studied  and  formulated  the 
formation,  growth  and  effect  of  a  secondary  vortex  in  the  near 
wake  of  a  circular  cylinder,  which  is  started  by  impulse  from 
rest  and  subsequently  moves  with  constant  velocity;  they  also 
studied  the  immediate  cause  leading  to  the  instability  of  the 
motion  of  symmetric  main  vortices  and  described  the  process  of 
the  foxmation  of  the  Karman  vortex.  Those  numerical  results 
are  given. 

In  this,  paper,  the  authors  suggest  that  under  the  effect 
of  a  secondary  vortex  while  the  symmetric  main  vortices  vary 
from  the  state  of  a  vorticity  supplement  into  free  vortices, 
the  final  state  has  no  further  vorticity  supplement.  The 
condition  that  leads  to  instability  of  the  motion  of  symmetric 
main  vortices  exists.  Then  undeT  the  action  of  a  small  distur¬ 
bance  a  fluctuating  motion  may  be  produced  and  further  develop¬ 
ment  leads  to  Karman  vortex  motion. 

In  the  present  work  the  flow  is  considered  to  be  laminar 
and  the  Reynolds  number  is  high.  The  results  given  in  the 
present  woTk  look  very  much  like  the  observed  features  in  flow 
field  visualization. 


I .  Foreword 


The  unsteady  motion  of  a  suddenly  moved  circular  cylinder  is  a  classic 
and  complex  problem.  With  a  high  Reynolds  number  and  laminar  flow  with  lengthen¬ 
ing  of  time,  these  complicated  flow  phenomena  include  unsteady  separation  of 
boundary  layers,  formation  and  growth  of  symmetric  vortices,  asymmetric  motion 
of  vortices,  formation  of  a  Karman  vortex,  vortex  street  motion,  and  others. 

Many  years'  theoretical  study  brought  about  very  successful  analysis  of  the  flow 
characteristics  within  a  very  short  time  after  start  of  motion  and  of  the  charac¬ 
teristics  of  vortex  street  motion  after  a  longer  time  of  motion;  fundamentally, 
the  study  solves  the  problem  of  an  unsteady  boundary  layer  under  the  condition  of 
potential  flow  solution.  However,  for  the  state  of  motion  in  the  intermediate 
stage,  that  is  the  stage  from  beginning  of  separation  to  the  formation  of  the 
Karman  vortex,  a  result  is  still  not  available  to  satisfy  the  theoretical  analysis 
except  a  few  purely  numerical  solutions  of  the  N-S  equations.  In  particular,  the 
following  three  problems  have  not  been  clarified:  (1)  the  formation  and  growth  of 
the  secondary  vortex  and  vortex  pair  as  well  as  their  effect  on  the  flow  and  main 
vortex;  (2)  direct  cause  and  condition  of  motion  instability  in  the  growth  process 
of  symmetric  main  vortices;  and  (3)  the  formation  process  of  the  Karman  vortex. 
Theoretically,  there  has  not  been  a  perfect  and  satisfactory  explanation  and  de¬ 
scription. 

Modern  experiments  utilize  a  delicate  flow  field  display  technique  to  conduct 
some  research  on  the  separation  process  of  secondary  and  main  vortices.  For  the 
main  work,  refer  to  papers  [1-4] .  Paper  [1]  was  the  first  observation  of  the 
existence  of  the  secondary  vortex,  and  it  pointed  out  that  Re-U»d/v>550  is  the 
condition  of  emergence  of  the  secondary  vortex.  In  paper  [1],  there  is  the  obser¬ 
vation  of  existence  of  a  pair  of  secondary  vortices  and  the  paper  derives  diagrams 
of  formation  and  breaking  up  .of  a  secondary  vortex.  The  experiment  in  paper  [2] 

clearly  indicates  the  formation  process  and  variation  of  the  main  vortex  near  the 
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wake  when  the  Re  number  is  10  to  2x10  with  formation  and  variation  of  a  second¬ 
ary  vortex.  It  was  observed  that  interruption  may  occur  in  the  main  vortex  layer 
at  a  certain  instant  in  the  growth  and  variation  process  of  secondary  and  main 
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vortices  when  Re>10;  no  interruption  occurs  in  the  main  vortex  layer  when  Re<10  . 
In  paper  [2] ,  it  was  pointed  out  that  the  state  of  the  secondary  vortex  is  closely 


rented  to  the  growth  and  destruction  of  main  vortex.  In  paper  [3]  (1979),  the 
experimental  results  when  the  Re  number  is  30-5000  were  reported.  When  Re>500, 
one  or  more  secondary  vortices  exist  and  have  been  observed  near  the  region  of 
separation  of  the  main  stream.  The  experiments  revealed  flow  variation  near  the 
vortex  in  the  formation  process  of  a  secondary  vortex.  These  experiments  enrich 
the  researchers'  knowledge  of  the  evolution  of  vortex  motion  behind  a  cylinder, 
which  carries  out  an  unsteady  motion;  this  is  helpful  in  regime  analysis  of  the 
flow  process.  However,  at  present  there  is  neither  a  corresponding  report  on  the 
theoretical  analysis,  nor  does  the  theoretical  calculation  result  match  the  experi¬ 
mental  result. 

Theoretically,  an  accurate  method  of  numerical  solution  was  used  to  calculate 
the  secondary  vortex;  the  main  work  was  done  in  papers  [5-8] .  In  paper  [5] ,  a 
numerical  solution  was  given  to  the  N-S  equation  under  a  high  Reynolds  number 
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(4x10  ) ;  the  formation  of  a  secondary  vortex  was  obtained  behind  the  cylinder. 

In  papers  [6-8] ,  calculations  were  conducted  for  a  medium  Reynolds  number  (Re 
about  500) ;  only  the  structure  of  a  single  secondary  vortex  was  derived  with  a 
blurred  diagram  of  slightly  bulging  streamlines  in  the  local  region.  There  is 
no  consistent  result  in  experimental  observations.  The  results  of  numerical  calcu¬ 
lation  only  matched  some  experimental  displays,  especially  in  formation  and  growth 
of  a  secondary  vortex.  Besides,  these  numberical  solutions  clearly  describe  the 
variation  process  of  the  flow  only  with  difficulty.  In  paper  [9] ,  a  model  of  a 
divergent  vortex  was  used  and  consideration  was  given  to  the  effect  on  flow  of 
separation  of  the  rear  shear  layer.  The  secondary  divergent  vortex  obtained  in 
the  calculation  enters  the  lower  stream  of  the  wake  with  the  structure  of  the 
vortex  blurred.  The  paper  did  not  conduct  a  theoretical  analysis  of  the  rear  shear 
layer  and  its  separation;  only  an  empirical  formula  of  the  separation  point  follow¬ 
ing  the  calculation  was  given.  Hence,  the  theoretical  study  conducted  achieved 
only  a  preliminary  result. 

This  paper  intends  to  use  a  method  of  combining  the  potential  flow  theory  of 
a  divergent  vortex  and  the  boundary  layer  theory.  Simultaneous  consideration  is 
given  to  the  model  of  a  divergent  vortex  separated  from  the  boundary  layer  and  the 
rear  shear  layer  in  order  to  study  and  calculate  the  formation,  growth  and  effect 
of  secondary  vortices  during  the  process  of  unsteady  motion  of  a  cylinder.  The 


direct  cause  and  conditions  are  also  studied  on  unsteady  notion  of  the  main  vor¬ 
tices,  and  the  process  of  formation  of  a  Karman  vortex  is  analyzed.  The  paper 
considers  the  situation  of  uniform-speed  straight-line  motion  after  a  cylinder 
suddently  moves  in  an  incompressible,  laminar  flow  with  high  Reynolds  number. 

II.  Theoretical  Analysis 

1.  Model  of  divergent  vortex 

In  paper  [10] ,  the  authors  described  in  detail  the  model  of  a  divergent  vor¬ 
tex.  Simultaneous  consideration  was  given  to  the  divergent-vortex  model  sepa¬ 
rated  from  the  boundary  layer  and  the  rear  shear  layer;  these  considerations  can 
be  generalized  as  follows:  the  boundary  layer  is  relatively  thin  with  a  high 
Reynolds  number;  the  stage  beginning  with  separation  to  the  formation  of  the  Karman 
vortex  is  the  initial  stage  of  the  entire  unsteady  motion;  in  this  stage,  a  region 
of  a  relatively  thin  shear  layer  of  circulation  flow  is  formed  gradually,  includ¬ 
ing  a  rear  stationary  point;  and  the  viscous  effect  is  mainly  concentrated  in  these 
thin  layers.  From  the  beginning  of  the  motion  to  the  instant  that  a  certain 
characteristic  occurs,  the  boundary  layer  bdgins  to  separate.  At  a  certain  time, 
under  certain  conditions  the  shear  layer  also  begins  to  separate.  The  separation 
points  vary  with  time  and  the  vortices  of  the  boundary  layer  and  rear  shear  layer 
move  to  the  lower-stream  wake  through  the  respective  separation  points.  The 
continuous  motion  process  of  this  vorticity  can  be  substituted  by  the  constant 
separation  and  motion  of  a  divergent,  nonviscous  potential  vortex  with  a  certain 
intensity  and  formation  at  certain  time  intervals.  However,  there  are  reverse 
directions  between  the  divergent  vortex  from  the  shear  layer  and  the  divergent 
vortex  formed  at  the  boundary  layer.  The  motion  of  the  divergent  vortex  is  deter¬ 
mined  by  the  local  potential-flow  field.  The  potential -flow  field  is  formed  by 
superposition  of  two  parts:  a  constantly  increased  divergent  vortex  and  the 
potential-flow  field  of  flow  winding  around  the  straight  uniform-flow  cylinder. 

With  time  variation,  flow  fields  corresponding  to  different  time  intervals  are  then 
obtained.  The  pressure  at  the  object  surface  and  velocity  distribution  are  new 
external-edge  conditions  of  the  boundary  layer  and  rear  shear  layer.  Under  these 
conditions,  separation  of  the  boundary  layer  and  rear  shear  layer  is  again  calcu¬ 
lated  to  determine  the  new  separation  point  and  divergent  vortex,  forming  a  new 
flow  field  of  the  entire  potential  flow.  By  repeated  calculation  in  such  a  way. 
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the  divergent-vortex  model  of  unsteady  initial  motion  with  high  Reynolds  number  is 
described  foT  a  cylinder  with  consideration  of  mutual  actions  among  boundary  layer, 
near-wake  flow  and  external  potential  flow. 


2.  Flow  field  of  nonviscous  potential  flow 

If  we  assume  that  W(Z)  is  the  compound  potential  of  potential -flow  field, 
and  r^  are  respectively  the  motion  positions  of  the  k-th  vortex  at  any  instant  of 
the  polar  coordinate  system,  then  in  paper  [10]  the  equations  are  as  follows: 
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In  the  equations,  Zn  and  Tn  are  respectively  the  position  and  intensity  of  the  n-th 

divergent  vortex.  N  is  the  total  number  of  all  divergent  vortices,  and  is  also  the 

summation  of  total  divergent  vortices  formed  with  separation  of  boundary  layer  and 

the  total  divergent  vortices  formed  by  separation  from  the  rear  shear  layer. 

Obviously,  N  varies  with  time.  Z_  is  also  a  function  of  time  t.  r  is  determined 

n  n 

by  vorticity  at  the  separation  point  of  the  boundary  layer  or  rear  shear  layer. 

The  equations  of  pressure  distribution  C^,  velocity  distribution  Vg  and  re¬ 
sistance  coefficient  at  the  object  surface  are,  respectively,  [10]: 
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3.  Separation  of  boundary  layer 

The  flow  closely  contacting  the  object  surface  should  satisfy  the  unsteady 
two-dimensional  incompressible  boundary  layer  equation  and  the  corresponding 
initial  boundary- value  conditions;  the  boundary  conditions  are  given  by  the  above- 
mentioned  distribution  of  velocities  of  potential  flow.  Hence,  the  internal 
layer,  external  layer,  and  near  wake  are  mutually  affected.  In  the  divergent 
vortex  model,  the  main  interest  is  the  variation  of  separation  point  with  time. 
Hence,  the  problem  can  be  regarded  as  the  coupling  solution  of  the  separation 
point  of  an  unsteady  boundary  layer.  Since  this  paper  emphasizes  the  problem  of 
the  secondary  vortex,  as  the  preliminary  description  of  this  problem  it  is  expected 
to  reduce  the  calculation  amounts  and  complexity  during  calculation.  Hence,  the 
variation  in  separation  point  of  the  boundary  layer  uses  an  empirical-relationship 
equation  in  paper  [10]  written  by  the  authors  and  for  the  time  being  the  coupling 
solution  is  not  given.  (This  empirical  equation  is  derived  from  an  accurate 
numerical  solution.  At  the  same  time,  the  calculated  numerical  results  of  the 
boundary- layer  separation  point  and  the  rear  shear- layer  separation  point  will 

be  given  in  another  paper.) 

9,  -  98.0  -  58.0«p[ -U968(<  -  0.65)]  (4) 

The  above  equation  indicates  the  variation  rule  with  time  of  a  separation  point 
for  actual  unsteady  flow  of  a  cylinder.  The  ini-  al  sepp  .ion  time  is  t*0.65 
and  0S«40*. 
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4.  Separation  of  reaT  shear  layer 

In  paper  [10]  a  detailed  analysis  was  given  of  the  flow  of  unsteady  motion  at 
its  initial  stage  without  considering  the  situation  of  the  secondary  vortex,  giving 
the  variation  rule  for  object- surface  flow  characteristics,  such  as  distributions 
of  object- surface  pressures  and  velocities.  On  this  basis,  this  paper  conducts 
another  step  of  regime  analysis  for  the  rear  shear  layer  formed  in  the  initial 
stage;  some  reasonable  assumptions  are  made. 

We  can  see  from  paper  [10] ,  on  the  formation  and  growth  of  the  symmetric  main 
vortices,  that  at  a  certain  instant  a  circulating  flow  region  appears  in  the  vicin¬ 
ity  of  the  object  surface  behind  the  cylinder.  With  lengthening  of  time,  the  range 
of  the  circulating  region  gradually  expands  and  includes  the  rear  stationary 
points.  At  that  time,  in  the  vicinity  of  the  object  surface  a  rear  shear  layer 
including  a  rear  stationary  point  is  formed.  From  the  rear  stationary  point,  the 
circulating  flow  along  the  object  surface  is  in  the  process  of  velocity  reduction 
and  pressure  increase  within  a  considerable  range.  This  paper  assumes  that  in  the 
initial  stage  of  this  unsteady  motion,  the  shear  layer  has  the  boundary  layer  char¬ 
acteristic  of  laminar  flow.  Under  the  persistent  (non-instantaneous)  action  of  a 
reverse  pressure  gradient,  boundary- layer  separation  occurs.  On  the  cylindrical 
object  surface  in  the  circulating  flow  region,  the  position  of  the  minimum  pres¬ 
sure  point  is  considered  to  separate  where  it  is  not  continuously  moving  upward 
along  the  cylindrical  surface.  At  that  time,  we  assume  that  the  variation  of  the 
rear  shear  layer  tends  toward  a  quasi -steady  state.  This  paper  uses  the  separa¬ 
tion-point  calculation  method  (on  the  rear  shear  layer)  for  the  boundary  layer  of 
Stratford  laminar  flow.  At  the  rear  separation  point,  the  following  relationship 
is  satisfied  [10] : 


*'-0-0, 


In  the  equation,  pmin^  is  the  minimum  pressure  value  on  the  cylindrical  surface 
in  the  circulating  flow  region;  9^  is  the  position  of  this  point  departing 


from  the  rear  stationary  point,  and  V*  is  the  velocity  of  external-edge  potential 
flow  of  the  corresponding  object  surface. 


I 

5.  Intensity  of  divergent  vortex 


The  intensity  of  the  divergent  vortex  should  equal  the  vorticity  of  the 

boundary  layer,  which  is  within  one  time  step  from  the  separation  point;  the 

intensity  can  be  indicated  by  the  velocity  of  potential  flow  of  the  object  surface, 
f  r 

Let  rn  and  rn  indicate,  respectively,  intensities  of  divergent  vortices  coming 
from  the  boundary  layer  and  rear  shear  layer;  and  Vqts  are,  respectively, 

velocities  of  external-edge  potential  flows  corresponding  to  the  separation  point 
of  the  boundary  layer  and  the  separation  point  of  the  shear  layer;  and  At  is 
the  length  of  the  time  step  during  calculation.  Then  [10] 
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4* 


III.  Brief  Description  of  Calculations 

Equations  (2) ,  (4) ,  (5)  and  (6)  are  simultaneous  equations  for  determining 
motions  and  intensities  of  N  divergent  vortices.  The  simultaneous  differential 
equations  are  numerically  solved  using  the  Runge-Kutta  method  of  varying  the 
number  of  equations.  In  calculations,  attention  should  be  paid  to  selecting  a 
unit  time  step.  If  a  time  step  is  too  large,  the  use  of  the  "divergent"  method 
to  simulate  vorticity  or  continuous  motion  of  vorticity  will  not  lead  to  accurate 
results.  If  the  time  step  is  too  small,  the  interval  of  divergent  vortices  will 
be  very  dense,  and  then  singularities  will  appear  while  using  the  theory  of 
potential  flow  to  calculate  the  induction  speed;  in  addition,  the  computer  time 
required  will  be  considerable.  The  paper  refers  to  paper  [10]  to  finally  select 
0.125  as  the  length  of  the  time  step.  For  the  time  of  initial  separation,  the 
paper  also  refers  to  paper  [10]  by  selecting  t*0.65.  The  newborn  vortex  is 
situated  at  9S*40°.  Later  on,  for  each  step  length  two  divergent  vortices  increase 
with  increase  of  six  equations.  (When  equation  (1)  is  used  to  calculate  the 
compound  potential,  the  contribution  of  the  display- image  vortex  is  included.) 


At  t=2.650  the  calculation  begins  on  the  separation  of  the  rear  shear  layer. 

Later  on,  for  each  step  length  four  divergent  vortices  increase  with  the  addi¬ 
tion  of  12  equations.  When  numrical  calculation  is  conducted  on  each  step  length 
of  equations  (2)  to  (6),  all  divergent  vortices  separated  from  the  boundary  layer 
and  all  divergent  vortices  separated  from  the  rear  shear  layer  should  be  taken 
into  account  and  numbered  consecutively.  After  calculation,  a  program  is  used  to 
distinguish  the  divergent  vortex  region  from  the  rear  shear  layer  already  existing 
at  that  instant  in  order  to  obtain  a  clear  image  of  the  secondary  vortex.  At  the 
separation  point  of  the  newborn  vortex  during  separation,  the  distance  (of  the 
vortex)  from  the  cylindrical  surface  should  not  be  greater  than  the  thickness  of 
the  boundary  layer  or  the  thickness  of  the  shear  layer.  For  calculation  conve¬ 
nience,  the  paper  selects  an  approximate  figure,  rn=1.05  [10],  The  paper  calcu¬ 
lates  from  separation  of  boundary  layer  to  tM.  In  the  final  calculation  of  step 
length,  there  are  74  divergent  vortices  and  222  simultaneous  equations,  including 
148  conventional  differential  equations.  In  order  to  save  computer  time  and  reduce 
internal  memory  storage  during  repeated  calculations,  the  calculation  results  of 
the  unsteady  process  are  stored  sector  by  sector,  and  the  compiled  program  can 
have  an  arbitrary  number  of  vortices  as  the  initial  condition  to  begin  the  calcu¬ 
lation.  By  using  equation  (3) ,  the  velocity  and  pressure  distributions  along 
the  object  surface  and  resistance  coefficient  can  be  obtained. 

IV.  Results 

As  revealed  by  analysis  and  theoretival  calculation,  in  the  unsteady  process 
of  a  cylinder  after  its  sudden  start  and  subsequently  maintaining  uniform  velocity, 
rear  shear  layers  are  formed  under  the  influence  of  induction  by  a  pair  of  main 
vortices  behind  the  cylinder.  Under  the  persistent  action  of  a  reverse  pressure 
gradient,  the  rear  shear  layers  are  separated  and  the  reverse-direction  secondary 
divergent  vortex  is  formed.  With  lengthening  of  time,  the  secondary  divergent 
vortex  forms  a  secondary  vortex  surface  within  a  limited  range;  the  vorticity  is 
gradually  accumulated  and  intensified.  The  formation  and  growth  of  the  secondary 
vortex  causes  apparent  variation  of  flow  in  the  vicinity,  and  apparently  affects 
the  motion  of  the  divergent  vortex  separated  from  the  boundary  layer  of  the  incom¬ 
ing  flow,  thus  affecting  the  motion  of  vorticity  separated  from  the  boundary  layer 
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and  the  accumulation  of  vorticity  of  symmetrical  main  vortices.  When  this  process 
develops  to  a  certain  degree,  the  surface  of  the  main  vortex  will  break  up.  At 
the  broken  terminal,  the  surface  of  the  main  vortex  again  rolls  up,  becoming  the 
vortex  surface  of  the  third  vortex.  Within  a  limited  space  range,  this  third 
vortex  and  the  secondary  vortex  form  a  secondary  vortex  pair.  At  another  broken 
terminal,  the  main  vortices  become  a  pair  of  free  vortices  almost  without  a 
vorticity  supplement.  In  paper  [11],  the  symmetrically  arranged  free  vortex  pair 
has  been  proved;  its  motion  is  unsteady.  Hence,  the  symmetrical  main  vortices 
vary  from  a  non-free  vortex  with  vorticity  supplement  to  free  vortices  without 
vorticity  supplement;  this  means  that  the  motion  of  symmetrical  vortices  has 
unsteady  conditions.  Due  to  the  effect  of  the  secondary  vortex,  the  main  vortices 
become  free  vortices;  this  is  the  direct  cause  of  the  symmetrical  main  vortices 
having  unsteady  motion.  At  that  time,  pulsation  occurs  easily  with  a  little 
disturbance  and  asymmetrical  motion  begins.  After  the  main  vortices  have  asym¬ 
metrical  motion,  the  secondary  vortex  will  break  up  and  develop  into  a  Karman 
vortex. 

Figure  (la  to  If)  shows  the  partial  calculation  results  in  the  paper.  Since 
in  the  near  wake  region,  within  the  time  duration  considered,  the  vortex  motion 
is  symmetrically  developed,  thus  in  the  figure  only  the  state  of  variation  of  the 
upper  half  of  the  plane  is  given.  Figure  la  shows  the  situation  of  symmetric 
vortices  in  the  near-the-wake  region  before  separation  of  the  rear  shear  layer. 
Within  the  time  duration  shown  in  Fig.  lb,  there  are  three  secondary  divergent 
vortices  and  the  secondary  vortex  surface  begins  to  be  formed.  Starting  from 
Fig.  lc,  the  secondary  vortex  is  constantly  intensified  with  more  and  more  effect 
on  the  main  vortices.  The  secondary  vortex  causes  an  apparent  change  of  vorticity 
motion  from  boundary  layer  to  lower  stream  and  accumulation  of  vorticity  of  the 
main  vortex.  At  the  time  instant  indicated  by  Fig.  lc,  the  surface  of  the  main 
vortex  breaks  up.  Later  on,  the  vortex  surface  of  the  broken  terminal  rolls  up 
and  forms  a  secondary  vortex  pair  [Fig.  If].  At  that  time,  almost  all  vorticity 
from  the  boundary  layer  is  accumulated  in  the  newly  formed  third  vortex.  Once 
the  main  vortices  become  free  vortices,  they  begin  entering  an  unsteady  state. 
These  calculation  results  describe,  with  an  image,  the  formation  and  growth  of 
a  secondary  vortex,  as  well  as  the  process  of  the  main  vortices  becoming  free 
vortices  due  to  the  effect  of  the  secondary  vortex.  The  formation  process  (at 
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the  initial  stage)  of  the  Karman  vortex  is  described  with  images.  In  paper  [2], 
the  vortex  motions  of  the  near-wake  region  of  a  cylinder  with  small  Reynolds 
number  (Re <520),  medium  Reynolds  number  (520<Re<1000)  and  high  Reynolds  number 
(Re >1000)  are  given.  The  photographs  (la-lg)  show  results  of  a  group  (among  them) 
of  high  Reynolds  number  (Re=1220) .  From  the  photographs,  the  growth  of  the 
secondary  vortex,  the  breaking  up  of  the  main  vortex  surface,  and  process  of  forma¬ 
tion  of  a  secondary  vortex  pair  (referring  to  photographs  l(c-d)  of  photo  plate  I) 
in  the  unsteady  motion  of  a  cylinder  can  be  clearly  seen.  Later,  beginning  with 
photograph  Id,  the  motion  of  the  vortex  pair  gradually  changes  into  an  asymmetrical 
situation.  Photographs  2  (a  to  f)  of  photo  plate  I  show  the  result  of  another 
group  of  experiments  (Re*1200)^ .  It  was  observed  that  the  growth  of  the  second¬ 
ary  vortex  and  variation  of  the  main  vortex  are  the  same  as  the  process  described 
in  paper  [2] .  The  results  of  analysis  and  calculation  in  the  paper  are  very 
similar  to  these  experimental  results  with  high  Reynolds  number.  Positions  of 
secondary  vortex  pairs  obtained  by  calculation  in  the  paper  are  situated,  respec¬ 
tively,  at  9^1°  and  51°.  The  result  is  quite  consistent  with  the  positions  (9*>40o 
and  50°)  of  secondary  vortex  paiTS  indicated  in  experiments. 


1^This  group  of  photographs  was  selected  from  "Display  Study  of  Vortices"  by  Matsui. 

Tatsuya  ,  Report  of  Japanese  Fluid  Dynamics  Delegation's  Visit  to  China, 

1979,  Beijing. 
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V.  Conclusion 

In  this  paper,  simultaneous  considerations  are  given  to  models  of  a  divergent 
vortex  with  separation  of  the  boundary  layer  and  shear  layer  with  analysis  and 
calculation  of  a  secondary  near-wake  vortex  separated  by  unsteady  motion  of  a 
cylinder  with  high  Reynolds  number.  The  paper  describes  the  formation  of  a 
secondary  vortex  and  the  effect  on  main  vortices,  as  well  as  the  process  of  change 
from  symmetrical  main  vortices  to  free  vortices.  The  paper  also  raises  the  condi¬ 
tions  of  unsteady  motion  of  symmetrical  vortices.  The  evolution  of  flow  in  the 
process  of  formation  of  a  Karman  vortex  is  described  in  detail.  The  regime  analy¬ 
sis  of  formation  of  the  Karman  vortex  is  provided.  The  formation  of  the  Karman 
vortex  nay  have  different  complex  details  with  different  Reynolds  numbers.  In 
this  paper,  consideration  is  given  to  the  flow  of  unsteady  motion  of  a  sector 
of  an  unsteady  process  from  a  sudden  start  to  the  appearance  of  asymmetrical 
motion  of  the  main  vortices;  the  paper  considers  the  conditions  of  high  Reynolds 
number  and  laminar  flow.  Further  study  is  required  in  theoretical  calculations 
from  the  sudden  appearance  of  asymmetrical  conditions  to  alternate  separation  of 
main  vortices  and  the  formation  of  a  Karman  vortex  row.  In  order  to  reduce  the 
calculations,  in  the  paper  the  variation  rule  of  separation  point  of  the  boundary 
layer  is  given  by  an  empirical  relationship  equation.  In  a  later  work,  the 
authors  obtained  the  numerical  calculation  results  for  the  solution  of  boundary 
layer  separation  and  separation  coupling  of  the  rear  shear  layer.  In  addition. 
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consideration  of  effect  of  a  viscous  vortex  eye  is  included.  Everything  concern¬ 
ing  the  formation,  growth  and  effect  of  the  secondary  vortex  and  the  variation 
from  main  vortex  to  free  vortex  is  consistent  with  the  descriptions  in  the  paper. 
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The  numerical  method  for  computing  the  viscous  compressible 
steady  laminar  flow  in  turbomachinery  [1]  is  applied  into  rota¬ 
tional  coordinate  system-  in  this  paper.  Fundamental  equations 
of  the  viscous  compressible  flow  in  the  rotor  are  introduced. 
Then,  the  numerical  solution  of  the  viscous  axisymmetrical  flow 
in  the  rotor  is  obtained  by  iterative  computation  based  on  the 
enthalpy  gradient  equation,  energy  equation  and  entropy  equation, 
expressed  along  arbitrary  curved  lines. 


I .  Introduction 


Since  Professor  Wu  Chunghua  proposed  a  flow  surface  theory  [2] ,  rapid 
development  was  made  in  the  calculation  method  of  three  dimensional  flow  in 
turbomachinery.  However,  in  these  calculations  the  fluid  is  considered  a  non- 
viscous  ideal  fluid;  the  mechanical  loss  caused  by  fluidity  of  the  fluid  can  be 
included  by  using  a  loss  factor  obtained  from  experiments.  Hence  the  calculation 
introduces  a  certain  error. 


Later,  Professor  Wu  further  proposed  a  viscous  flow  theory  [3]  of  the  turbo¬ 
machinery  flow  channel;  a  detailed  analysis  was  conducted  on  the  physical  signif¬ 
icance  of  viscous  items.  These  studies  laid  a  foundation  for  the  numerical  calcu¬ 
lation  of  viscous  flow  in  turbomachinery.  In  recent  years,  Patankar  and  Spalding 
[4]  used  the  progressive  integration  method  to  calculate  viscous  incompressible 
flow  in  equal-cross-section  straight  pipes.  Pratap  and  Majundar  [5]  used  the 
same  method,  and  calculated  viscous  incompressible  flow  of  rectangular  cross- 
section  straight  pipes  with  flow  rotating  around  the  axis  perpendicular  to  the 
direction  of  the  main  stream.  However,  the  calculation  of  viscous  compressible 
flow  problems  in  the  rotor  flow  channel  of  turbomachinery  did  not  lead  to  any 
solution  up  to  now. 

In  the  first  part  of  paper  [1] ,  the  author  proposed  a  numerical  method  of 
calculating  the  viscous  compressible  steady  laminar  flow  problems  of  an  internal 
flow  field.  In  addition,  numerical  calculations  were  conducted  on  viscous  axisym- 
metric  flow*  of  a  turbomachinery  stator.  The  paper  extends  this  method  into  a 
rotating  coordinate  system,  and  derives  fundamental  equations  satisfying  viscous 
compressible  flow  in  a  rotating  coordinate  system  with  rotation  around  a  fixed 
axis  at  constant  angular  velocity.  Iterative  calculation  using  these  equations 
can  lead  to  numerical  solution  of  viscous  compressible  axisymmetric  flow  in  the 
rotor  channel  of  turbomachinery.  This  can  be  approximately  considered  the  numeri¬ 
cal  solution  on  the  intermediate  flow  surface  S2-  On  this  basis,  similar  methods 
are  used  to  conduct  numerical  calculation  in  the  flow  surface  Sj.  Then  the  numeri¬ 
cal  solution  of  quasi  three  dimensional  flow  in  the  rotor  can  be  obtained. 

II.  Fundamental  Equations  in  Rotating  Coordinate  System 

According  to  the  characteristics  of  streams  in  turbomachinery,  a  rotating 
coordinate  system  with  rotation  along  the  rotor  can  be  established.  If  it  is 
assumed  that  the  rotating  angular  velocity  of  the  rotor  is  constant  and  the  stream 
in  the  rotor  is  an  ideal  gas,  which  satisfies  Newton's  viscosity  law  and  Fourier's 
law  of  thermal  conduction,  then  as  given  in  the  appendix  using  the  (flow  parame¬ 
ters)  transformation  equation  between  a  stationary  coordinate  system  and  rotating 
coordinate  system,  it  is  possible  to  process  the  fundamental  equations  of  viscous 

*Re^er  to  paper  [1]  for  mechanical  significance  of  axisymmetric  flow. 
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compressible  steady  flow  J.n  a  stationary  coordinate  system  in  order  to  derive 
fundamental  equations  satisfying  the  viscous  compressible  relative  steady  flow  in 
the  rotating  coordinate  system: 

V  •  (pw)  -  0  (1) 

(10  •  V)u»  +  2a»  x  nr  +  X  (<*  X  A)  —  F—~Vp  +  —  V  •  [r]  (2) 

p  P 

m-  vh  —  Q  —  —  V  •  *+  —  v  Vp  +  i-  [r]  •  Vw  (3) 

P  P ,  P 

p-RpT  ,  .  ,  (4) 

k-c,T  (5) 

In  the  equations,  w  is  the  relative  flow  velocity  of  the  stream;  p,  p,  T  and  h 
are  (respectively)  the  pressure,  density,  temperature  and  enthalpy;  R  is  the  gas 
constant;  c^  is  the  specific  heat  at  constant  pressure,  F  is  the  mass  force;  Q 
is  the  amount  of  heat  absorption;  to  is  the  rotating  angular  velocity;  R  is  the 
position  vector  of  the  stream;  [t]  is  the  viscous  stress  tensor;  and  q  is  the 
amount  of  heat  conduction.  These  parameters  can  be  expressed  as 

[r]  - p\vw  +  (V»)T -  J-  v  • «(/}}  («) 

(7) 

T 

In  the  equations,  (Vw)  is  the  transposed  tensor  of  tensor  (vw);  [I]  is  the  unit 
tensor;  and  y  and  <  are,  respectively,  the  viscous  coefficient  and  heat  conduction 
coefficient  of  the  fluid.  When  the  temperature  varies  within  a  narrow  range, 

Cp,  y  and  <  can  be  approximately  considered  as  constants. 

The  unknown  functions  (w,  p,  p,  T  and  h)  are  included  in  equations  (1)  to 
(5);  therefore  the  equations  aTe  closed.  Before  using  this  method  to  conduct 
numerical  calculation,  it  is  necessary  to  first  conduct  transformation  of  the 
equations  by  using  the  similar  methods  in  paper  [1].  Also,  the  transformed 
fundamental  equations  in  paper  [1]  can  be  applied  by  using  the  transformed  rela¬ 
tionship  equations  of  flow  parameters  in  two  coordinate  systems  in  the  appendix 
in  ordeT  to  derive  the  following  transformation  equations  in  the  rotating  coordi¬ 
nate  system: 


Flow  equation 


fiW-JA  —  G9 


C») 


In  the  equation,  A  is  the  area  of  any  cross  section  in  the  flow  channel,  and  Gq 
is  the  given  flow. 


The  enthalpy  gradient  equation 

VA  —  TV/  +  F+  V  •  [r]  —  («*  •  v)ur  —  2m  X  m  —  m  X  (m  x  R)  (9) 

In  the  equation,  s  is  the  entropy  of  the  fluid  with  the  equation  of  entropy  as 
follows:  * 

'  “  R^~-  (10) 

In  the  equation,  K  is  the  adiabatic  index  of  the  fluid. 


The  energy  equation 

^f-,7(p+ » v‘T+7',-ar]-->}  a» 

In  the  equation,  l  is  the  spatial  streamline;  d/dt  is  the  direction  derivative 
along  the  streamline  t;  I  is  the  stagnant  transposed  enthalpy  and  its  equation  is 

1  “  *  +  ™  +  *2  ”  rmu  “  **  “  02) 

In  the  equation,  v  is  the  absolute  flow  velocity  of  the  stream  and  hQ  is  the 
stationary  enthalpy  of  the  fluid. 

The  entropy  equation 

-  -V  {  0  +  4 ^  +  T  [r] '  03) 

dl  wT  \  P  P  • 

Then,  the  transformed  equations  are  composed  of  equations  (4),  (5),  and  (8) 
to  (13).  In  these  equations,  equation  (9)  along  the  streamline  direction  can  be 
used  to  derive  energy  equation  (11) .  Therefore,  only  two  directions  can  be 
applied  in  equation  (9).  In  the  equations,  the  unknown  functions  included  are 
w,  P»  P*  T,  h,  s  and  I;  therefore,  the  equations  are  closed  and  can  be  used  to 
conduct  numerical  calculation. 


III.  The  Boundary  Conditions  of  Equation  Forms  in  Turbomachinery 

In  turbomachinery,  generally  the  mass  force  F  and  heat  absorption  Q  are  not 
considered.  However,  in  the  axial  symmetric  flow  discussed  in  the  paper,  the 
partial  derivative  to  6  of  all  flow  parameters  is  0.  Then,  the  above-mentioned 
equations  can  be  written  in  the  foim  of  a  cylindrical-surface  coordinate  system. 

The  enthalpy  gradient  equation 

For  any  spatial  curve  q  in  a  flow  field,  the  directional  derivative  of 
enthalpy  h  along  the  direction  of  curve  q  is 

{"7  V  “  (■*  "  *)*■ -2mx  m  —  mx  (»x  it)}  (14) 

In  the  equation,  q°  is  the  unit  vector  along  the  direction  of  curve  q.  Here,  q 
can  be  assumed  as  three  directions  not  on  the  same  plane,  thus  deriving  an  equation 
satisfying  the  three  directions.  First,  a  direction  of  streamline  1  is  applied; 
the  equation  along  this  direction  was  used  to  derive  energy  equation  (11).  The 
6  direction  is  taken  as  the  second  direction;  for  the  axisymmetric  flow  discussed 
in  the  paper,  this  equation  does  not  need  to  be  discussed.  The  third  direction 
is  the  direction  of  the  station  line  [1]  approximately  perpendicular  to  the 
direction  of  either  stream  l  or  8.  Then,  the  enthalpy  gradient  equation  along 
the  station  line  can  be  written  as 
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In  the  equation,  n  indicates  the  projection  of  the  station  line  on  the  meridian 
plane;  this  is  the  meridian  station  line,  d/dn  indicates  the  total  derivative  of 
flow  parameters  varying  along  the  spatial  station  line  to  the  meridian  station 
line  n.  d/dm  indicates  the  total  derivative  of  the  meridian  streamline  m  of  flow 
parameters  varying  along  the  spatial  streamline,  w  ,  w ,  and  w  are  components  of 
the  relative  velocity  w  in  the  cylindrical -surface  coordinate  system  (r,  0  and  z) ; 
wffl  is  the  projection  of  w  on  the  meridian  plane.  Among  them,  the  following 
equations  are  satisfied. 


m/#  —  wtinfl,  mm  —  mcotf,  I 


oo 


In  the  equation,  a  is  the  meridian  flow  angle  and  6  is  the  circumferential-direc¬ 
tion  flow  angle  (Refer  to  Fig.  1) . 


Fig.  1. 


For  simplicity  in  calculation,  generally  the  meTidian  station  line  n  can  be 
taken  as  a  straight  line;  this  is  the  quasi-perpendicularly  intersecting  line  as 
it  is  usually  called.  For  axial-flow  turbomachinery,  the  meridian  station  line  n 
can  be  more  simply  applied  as  the  diametral-direction  station  line.  At  that  time, 
(dz/dn)=0  in  equation  (15),  which  can  be  simplified  into 
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In  the  equation,  d/dr  4s  the  total  derivative  of  diametral  direction  r  of 
flow  parameters  varying  along  the  spatial  station  line. 

The  energy  equation 
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The  entropy  equation 


is  -( 

dm 


+ 2  (*y + & -f-r  *  (i?y + & + ^y 


(18) 


All  the  three  above-mentioned  equations  are  written  in  derivative  form  along 
the  non-perpendicular ly- intersected  curve  directions  m  and  n.  In  explicit  itera¬ 
tive  calculation,  these  equations  are  simple  differential  equations.  Numerical 
calculation  can  be  conducted  through  numerical  integration. 

For  numerical  calculation  of  viscous  compressible  axisymmetric  flow  in  a 
turbomachinery  rotor,  it  is  necessary  to  apply  the  following  boundary  conditions: 
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1.  On  the  internal  and  external  wall  surfaces  in  the  rotor  flow  channel,  the 
relative  velocity  should  be  0;  in  other  words,  wwan=0-  The  temperature  on  the 
wall  surface  can  be  given  in  two  ways:  (1)  with  the  temperature  distribution  on 
the  wall  surface  given,  and  (2)  treated  as  an  adiabatic  wall  surface,  that  is 
the  normal -direct ion  derivative  of  the  temperature  on  the  wall  surface  is  0. 
According  to  the  relationship  between  temperature  and  entropy  satisfying 
equation  (5),  the  boundary  conditions  of  enthalpy  on  the  wall  surface  can  be 
correspondingly  written. 

2.  When  calculating  the  rotor  flow,  in  order  to  eliminate  the  error  of  arbitrarily 
selected  boundary  conditions  at  inlet  and  outlet,  it  is  necessary  to  extend  a 
distance  toward  the  upper  and  lower  streams.  As  there  is  no  function  of  blade 
force  in  the  extension  region,  in  equation  (14)  q  can  be  taken  as  0,  and  from 
(ah/a0)=O,  we  derive 

4Lr.(y  +  £«0J  — f  +  i  Qy  4  __  1  /•,  q\ 

dm  pwm  L  dr1  r  dr  d*  U ' 

From  the  above  equation,  the  distribution  of  wQ  in  the  extension  Tegion  can  be 
calculated.  On  the  wall  surface  is  the  extension  region,  wq  can  be  given  in 
two  ways:  (1)  if  the  wall  surface  of  the  extension  region  rotates  along  with  the 
rotor,  it  can  be  derived  that  we|wau®0*  (2)  If  the  wall  surface  of  the  extension 
region  remains  stationary,  it  can  be  given  that  wq  j  Wa.ll=-ru> * 

IV.  Numerical  Calculation  Method  and  Calculation  Examples 

This  paper  adopts  the  iterative  substitution  method  of  the  streamline  for 
numerical  calculation;  the  calculation  steps  are  basically  the  same  as  the  calcu¬ 
lation  steps  mentioned  in  paper  [1];  some  slight  differences  exist  in  some  local 
calculations.  The  calculation  steps  of  axisymmetric  flow  of  a  rotor  are  briefly 
described  as  follows: 

1.  On  the  intermediate  flow  surface  S2  (including  the  inner  region  of  the  rotor 
and  the  extension  region  toward  the  upper  and  lower  streams) ,  several  station 
lines  q  are  selected.  According  to  certain  distribution  principle  of  flow,  the 
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assumed  streamline  l  is  marked.  The  streamlines  and  station  lines  compose  the 
calculation  network;  on  network  points,  the  distribution  situation  of  v  ,  p  and  T 
are  assumed. 

2.  The  calculation  network  is  projected  onto  the  meridian  plane  to  obtain  a 
meridian  network  composed  of  meridian  streamline  m  and  meridian  station  line  n. 

By  using  numerical  differentiation,  the  meridian  flow  angle  a  and  meridian  station 

line  angle  »  — acts  —  ,  are  obtained. 

dr  \  along  n 

dd 

3.  By  using  an  arc  shape  of  the  blade,  we  calculate  - 

dm 

For  boundary  condition  (19),  iterative  calculation  is  conducted  to  obtain  w  in 

0 

the  extension  region  of  the  upper  and  lower  streams.  To  reverse  the  problem, 
w  can  be  calculated  according  to  the  selected  flow  type.  Then,  we  obtain  the 
circumference-direction  flow  angle  £  —  Mcig-21, 

»m 

4.  By  using  direction  derivative  equations  (23)  and  (24)  in  paper  [1],  the  partial 
derivatives  of  the  velocity  component  and  temperatures  of  r  and  z  can  be 
calculate:. 

5.  <fci  each  streamline  (except  the  wall  surface),  we  conduct  numerical  integration 
for  energy  equation  (17)  and  entropy  equation  (18)  to  obtain  distributions  of 
stagnent,  transposed  enthalpy  and  entropy  in  the  flow  field. 

6.  On  each  station  line,  we  conduct  numerical  integration  of  the  enthalpy  gradient 
equation  (IS).  In  addition,  by  using  the  boundary  condition  of  the  wall  surface, 
the  distribution  of  enthalpy  h  can  be  obtained.  Then  by  using  equations  (12),  (5), 
(10)  and  (4),  the  values  of  w,  T,  p  and  p  can  be  obtained.  Then  by  using  the 
flow  equation  (8),  the  numerical  integration  constant  can  be  determined. 

7.  According  to  the  pre-determined  flow  distribution  principle,  new  streamline 
positions  can  be  obtained  by  reverse  calculation.  Then,  we  return  to  step  2 

to  conduct  the  iterative  calculation  of  the  next  cycle  until  the  required  accuracy 
is  satisfied.  At  that  time,  the  calculation  result  is  the  required  numerical 
solution. 


along  streamline 
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The  characteristics  of  this  method  are  simple  and  brief  in  principle,  with 
fast  calculation  speed,  and  not  too  much  computer  internal  storage  required. 
Hence,  this  method  can  be  extensively  applied  to  engineering  design. 


Appendix:  Transformation  Equations  Between  the  Stationary  and  Rotating  Coordinate 
Systems  of  Plow  Parameters  of  a  Fluid 


For  a  rotating  coordinate  system  established  in  the  papeT,  the  following 
relationships  exist  between  the  stationary  and  rotating  coordinate  systems  of 
flow  parameters. 


1.  The  parameters  of  the  flow  state  do  not  rely  on  the  selected  coordinate  system. 
The  total  derivative  of  time  and  the  partial  derivatives  of  various  coordinates 
also  do  not  rely  on  the  selected  coordinate  system,  that  is 


it  £L. 

4t  4l 

en  at  «rt  at  *7 

«r“ar»  «i"  h 


(Appendix  1) 
(Appendix  2) 


In  the  equations,  f  is  the  parameter  of  a  certain  state.  The  apostrophe  at  the 
upper  right  corner  of  the  differential  symbol  indicates  that  the  derivation  process 
is  conducted  in  the  rotating  coordinate  system.  Otherwise,  the  derivation  process 
is  conducted  in  the  stationary  coordinate  system.  In  other  parts  of  the  paper, 
as  all  discussions  are  conducted  in  the  rotating  coordinate  system,  this  apostrophe 
is  omitted. 


2.  In  a  flow  field,  an  arbitrary  vector  o  (for  example,  the  position  vector  R 
of  the  mass  point  of  the  fluid)  does  not  depend  on  the  selected  coordinate  system. 
The  total  derivatives  of  time  are  different  in  two  types  of  coordinate  systems, 
that  is 


4m 

4$ 


4‘m 

"  ~7T  *  * 


(Appendix  3) 


If  a  is  substituted  by  R,  the  above  equation  can  be  written  as 

"L  m .£*  +  w  x  A, 


that  is,  •  »a’+«x*  (Appendix  4) 

Or,  it  is  written  in  component  form 

*,  -  *„  r,  ■»  »,  +  rm,  r,  -  m,  (Appendix  4  ' ) 
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3.  By  using  equations  (appendix  3)  and  (appendix  4),  the  relationship  between  the 
absolute  acceleration  dv/dt  and  the  relative  acceleration  d'w/dt  in  the  rotating 
coordinate  system  can  be  derived: 


(•**) 


(Appendix  5) 


4.  According  to  the  above-mentioned  equations,  the  gradients,  divergences,  degxees 
of  rotation,  and  Laplace  calculation  symbols  of  an  arbitrary  scalar-quantity  field 
$  or  vector  field  a  in  the  flow  field  can  be  derived.  (None  of  these  parameters 
rely  on  the  selected  coordinate  system.)  That  is, 

-  v>,  tx«- (Appendix  6) 


S.  By  using  the  above-mentioned  results,  the  following  results  can  also  be  derived: 


▼  •  W  -  ▼'.•  I*']*  1*1  • w 


(Appendix  7) 


In  the  equations,  refer  to  paper  [1]  for  expression  [T];  refer  to  equation  (6)  in 
the  paper  for  expression  [t • J . 


LITERATURE 


1.  Jiang  Jinliang,  Numerical  Calculation  of  Viscous  Compressible  Flow  in  Turbo¬ 
machinery  Flow  Channel.  (I)  Calculation  of  Axisymmetric  Flow  in  a  Stator,  LIXUE 
XUEBAO  [ACTA  MECHANICA  SINICA] ,  2  (1981) . 

2.  Wu  Chunghua,  A  General  Theory  of  Three-Dimensional  Flow  in  Subsonic  and 
Supersonic  Turbomachines  of  Axial,  Radial  and  Mixed-Flow  Types,  NACA  TN  2604, 

1952. 

3.  Wu  Chunghua,  Fundamental  Equations  of  Aerodynamic  Thermodynamics  in  Stationary 
and  Dynamic  Coordinate  Systems  —  Action  of  Viscous  Force  and  Physical  Significance 
of  Viscous  Items,  JIXIE  GONGCHENG  XUEBAO  [JOURNAL  OF  MECHANICAL  ENGINEERING],  13, 

4  (1965). 

4.  Patankar,  S.  V.,  and  Spalding,  D.  B.,  A  Calculation  Procedure  for  Heat,  Mass 
and  Momentum  Transfer  in  Three-Dimensional  Parabolic  Flows,  INTERNATIONAL  JOURNAL 
OF  HEAT  AND  MASS  TRANSFER,  15  (1972) . 

5.  Pratap,  V.  S.,  Majunder,  A.  K.,  and  Spalding,  D.  B.,  Numerical  Computation  of 
Flow  in  Rotating  Ducts,  T.  ASME.  SI.,  99,  1  (1977). 


27 


6.  Jiang  Jinliang,  Streamline  Iterative  Substitution  Method  of  Calculating 
Two-Dimensional  Viscous  Flow,  SHANGHAI  LIXUE  [SHANGHAI  MECHANICS],  2,  1  (1981) 

7.  Koch,  C.  C. ,  Bilwakesh,  K.  R. ,  and  Doyle,  V.  L. ,  Evaluation  of  Range  and 
Distortion  Tolerance  for  High  Mach  Number  Transonic  Free  Stages,  NASA  CR- 72806 
I,  II  (1971). 


PARTIAL  DIFFERENTIAL  EQUATIONS  SATISFIED  BY  STREAM  FUNCTIONS  IN  THREE  DIMENSIONAL 
TURBOMACHINERY  FLOW  AND  THE  FINITE  ELEMENT  SOLUTION 


Li  Kaitai ,  and  Huang  Aixiang 
Jiao tong  University,  Xi'an 
Submitted  September  1980 

By  applying  a  semi-geodesic  coordinate  system  and  tensoT 
analysis,  this  paper  deduces  partial  differential  equations 
satisfied  by  stream  functions  on  any  stream  surface  of  turbo¬ 
machinery  flow.  The  paper  has  also  defined  a  weak  problem 
associated  with  the  boundary  value  problem  for  a  quasilinear 
elliptic  partial  differential  equation,  and  has  provided  an 
algorithm  of  a  finite  element  solution.  Two  examples  are 
calculated.  The  establishment  of  this  partial  differential 
equation  has  brought  about  a  unified  method  to  find  the  flow 
solution  on  any  stream  surface.  Therefore  a  unified  program 
can  be  established  to  realize  alternative  iteration  between 
Sj  and  Sj  on  the  basis  of  Wu's  theory. 

I.  Flow  Surface  and  Its  Measured  Tensor,  as  Well  as  Single-Parameter  Stream- 
Surface  Family  Covering  Three-Dimensional  Flow  Field 

We  assume  any  coordinate  system  (x^)  in  the  three-dimensional  Euclidean  space 
R3.  If  the  parameter  equation  of  a  two-dimensional  curved  surface  S  is 

.-1,2,3 


In  the  equation,  u  (a»l,  2)  is  the  Gaussian  coordinate  system  of  S.  Then  the 

measured  tensor  g. .  of  R'  below  (x^)  and  the  measured  tensor  a  .  of  curved  surface 
ij  08 

S  have  the  following  relationship: 

-  gtfiJ/dtf  •  da*/dm>  (U) 

The  reverse  variation  component  is  —  «* 

That  is, 

*“  *tJ»  (1.2) 

* 

In  the  equations,  a*| a^ |  is  the  matrix  of  the  measured  tensor. 

For  a  three-dimensional  flow  field  inside  the  turbomachinery,  we  assume  the 

12  3 

cylindrical  coordinate  system  y  »r,  y  *$,  and  y  »z  rotating  with  the  turbine.  At 
that  time,  the  measured  tensor  is 


ltn  1  • 


10  0 

0  r*  0 

0  0  1 


(U) 


If  there  is  a  flow  surface  of  internal  flow  in  a  turbomachinery  flow  channel, 
the  parameter  equation  is 

y'-y'O-)  (1.4) 

Then  the  measured  tensor  of  the  flow  surface  is 


W  dm-  dm?  0*  dm?  du'  dm*  dmm  du> 
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The  flow  surface  can  be  considered  formed  with  a  certain  rule.  Through  an 
arbitrary  point,  the  streamline 

™  (y* "  ri>  y»  —  <p»»  y»  —  *«) 

is  the  integration  curve  of  the  following  initial-value  problem 


it 

r'U-rt 


(1.0 


In  the  equation,  s  is  the  arc  length  of  the  streamline.  However 


f  -  WJW,  f  -  WjW/r,  f  -  WJW 

W  is  the  relative  flow  velocity-of  the  stream;  W«|W|,  W  .  W . ,  and  W  are  physical 

components  of  W  in  the  relative  cylindrical  coordinate  system.  The  solution  of 

equation  (1.6)  can  be  expressed  as 

yl-yi(s;  Yo)  (1.7) 

We  assume  that  Y_  varies  on  the  certain  curve  t  . 

o  o 

V.-Y/r)  r,  <  r  <  r.  (M> 

Then  equation  (1.7)  can  be  expressed  as 

>•-/(<;  y •(*■)) -so,  o  <l-9> 

It  is  easy  to  prove  that  in  case  no  characteristic  direction  exists  on  lQ,  the 
order  number  of  the  Jacobi  matrix  of  equation  (1.9) 


/ 


'*£*£*& 

At  A*  At 

At*  Ag4  Ag* 


(1.10) 


is  2.  Hence,  the  streamlines  passing  through  any  point  of  i  compose  a  curved 

1  °  2 

surface  S;  this  is  the  curved  surface.  In  addition,  s*u  and  T*u  can  be  consid¬ 
ered  the  Gaussian  coordinate  system  of  the  flow  surface.  Equation  (1.9)  is 

the  parameter  equation  of  the  flow  surface  with  the  measured  tensor  a  .  After 

ap 

substituting  equations  (1.9)  and  (1.6)  into  equation  (1.5),  we  derive 


•u  -  tiff 


gif  —  gt*/V4 


Af'  Og!  ,  ,  a  ri  A  »f 

6t~ 


(i.ii) 


In  the  equations,  the  matrix  — i 

equation  (1.6)  - 

"  fi  ('»  y)*i  ) 


is  the  solution  of  the  variation 


fi  ■*  Qf/dy*  and  5  *  are  the  Kronecker  signs . 

When  lQ  moves  with  certain  forms:  /» —  4(/),  Only  without 

appearance  of  the  characteristic  direction  at  an  arbitrary  position  t €  /j 
in  the  motion  process,  there  is  no  emergence  of  a  characteristic  direction  on  2,  . 
Then  the  integration  flow  passing  through  2q(1 .6)  forms  a  single-parameter  flow 
surface  family  S(t),  whose  Riemann  structure  is  given  in  equation  (1.11).  In 
addition,  (<»  *tl)  covers  the  entire  three-dimensional  flow 

field. 

For  example,  for  inlet  flows  without  rotation,  a  series  of  arc  sectors  from 
blade  root  to  tip  (at  the  turbine  inlet)  can  be  selected  as  <LQ 

- 0:»> 


At  that  time  f  but  •d_«,  m»_t,  •&- 0.  Equation  (1.11)  can 

0r  Or  Or 

be  expressed  as 

-  (1.14) 

while  taking 

KRt-Rj  9<i<  l 

Then  iQ(t)  moves  toward  the  blade  height  direction.  Here,  R1  and  R2  are  radii  of 
the  rotor  cover  and  rotor  at  the  inlet.  When  t([0,l]  varies,  S(t)  covers  the 
three-dimensional  flow  field  of  the  entire  flow  channel. 


—  come 

rf"f«  9i 

yi  conat 
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Key:  (1)  Extension  flow  channel  toward 
the  lower  stream;  (2)  Extension  flow 
channel  toward  the  upper  stream. 


II.  Differential  Equations  of  Flow  Functions  on  Arbitrary  Flow  Surfaces 

We  assume  that  a  turbine  rotates  with  angular  velocity  u>;  W  is  the  relative 

2 

flow  velocity  of  the  fluid;  T  is  the  fluid  temperature;  s  is  entropy;  I»i+(W  /2)- 
2  2 

U  r  /2)  is  the  transposed  enthalpy;  and  i  is  the  density  enthalpy  of  the  fluid. 
Then  the  continuous  equation  and  momentum  equation  of  internal  flow  of  the 
turbomachinery  are 

v  •  OWO  -  o  -  0.1) 

wy-imx  w— vi -tvs  (lay 

In  the  equations,  p  is  the  fluid  density. 

We  introduce  the  matrix  tensor 


V7 

*  • 

o 


l/Vi 

-wT 

o 


where  i,  j  and  k  are  even  row  of  (1,2,3) 
where  i,  j  and  k  are  odd  row  of  (1,2,3) 
for  other  situations 

(2.3) 


In  the  matrix,  g"|g.,  is  the  matrix  of  the  measured  tensor.  e 


ijk 


and  are 


asymmetric  to  any  pair  of  indexes.  In  addition,  there  is  the  accompanying  tensor: 


•irt  “  tlmgfrtot* 


(2.4) 


The  degree  of  rotation  and  vector  product  can  be  expressed  as 


VXW- -  7,,W!  1 

«  X  I 

\ 

In  the  equations,  {e. }  is  the  covariant  fundamental  vector  of  (xl)  coordinate 

i  1 

system;  and  w  are,  respectively,  the  covariant  component  and  reverse  variant 
component  of  W;  and  7^  is  the  covariant  derivative  of  the  tensor.  Therefore, 
equation  (2.2)  can  be  expressed  as 

—  -  ^(V,/  —  r7*s)  (2.«) 

2 

ik 

g  is  the  reverse  variant  component  of  the  measured  tensor. 


First,  the  semi-geodesic  coordinate  system  (x1)  of  the  stream  surface  S  is 
considered  as  the  fundamental  curved  surface;  then  the  following  relationships 
obtain  between  the  measured  tensor  g„  of  R3  below  (xi)  and  the  measured  tensor 
aa@  stream  surface  S: 

—  — 0*  fw— 1*  g  — •  f  . 

|W  mm  gtm  a  “  0»  |*™  | 

Under  the  semi-geodesic  coordinate  system, 

*>5|,  —  —  0  (2.8) 

We  substitute  the  equations  (2.7)  and  (2.8)  into  equation  (2.6),  and  we  obtain 


i-  —  V,»»)  - 

-  «'•(?./  -  TVs) 


(2.9) 


By  utilizing 
obtain 


^  eijk 


concerning  their  index  asymmetric  characteristic,  we 


y  «*^(Vj«*'i  —  Vtwx) 

“  (aT‘w*  —  —  Vtm>i) 


From  equation  (1.2)  and  *a  « equation  (2.9)  can  be  expressed  as 

w»I(V,«*  -  +  2  \/T  *»]  -  *«‘*(V./  -  7V.i)  1 

—  V,*,)  +  2  —  —s**(VJ—TVs)  ' 
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or 


+  2\/  « <*3  “  / 

f  ««“(?./  -  TV;)/*,  when  W2^0 
/  ”  !-««*(?./  -  TVj)/wt  when  w,^0 


(2.10) 


(2.11) 


Under  the  semi-geodesic  coordinate  system,  the  normal  direction  of  the 
stream  surface  is  n,  which  is  consistent  with  e. : 


*  —  •*  —  •*,  *••.  —  <) 

Then,  from  equation  (2.5)  we  obtain 

it  *  *ot Iff  — 

Substituting  into  equation  (2.10),  we  obtain 

i»  •  roc  IV  +  2 j*  •  —  //v"7 

t  -  « 

This  is  the  equation  of  the  degree  of  rotation  on  the  stream  surface. 
Equation  (2.1)  can  be  expressed  as 


(Vw)  +  ~  (v'  •  * 


By  noting  equation  (2.8),  on  the  stream  surface  there  is 


(2.12) 


(2.13) 


(2.14) 


— +  p  -/•  -  o 


(2.15) 


We  assume  b(x  )  satisfies  an  one-order  partial  differential  equation: 

flu*  I  _  dm* 

I.  dx*  >  (2.16) 

Then,  we  obtain  a  continuous  equation  (on  the  stream  surface) 


(2.17) 


In  the  equation,  b(x°)  is  an  arbitrary  integration  of  equation  (2.16).  It  can 
be  proved  that  the  mechanical  significance  is  the  relative  thickness  of  the  stream 
slice. 

By  introducing  flow  function  ^(x01)  satisfying 
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dj/dx*  —  —  bp^mm?  |  C2.18) 

d<t>/dx*  ^  bp  */  *  wl 

/ 

If  exists,  then  wa  determined  by  equation  (2.18)  can  automatically  satisfy  the 
continuity  equation  (2.17). 

By  substituting  equation  (2.18)  into  equation  (2.10),  we  obtain  the  equation 
(of  Hx01)) 


In  the  equation, 
curved  surface. 


ftr*  V  bp  d*J  ■ 


(2.19) 


p  mm  2  >/"*"“*  —  1-  When  the  stream  surface  is  the  real  Riemann 


ag 

a  is  positively  symmetric. 


The  boundary  I’  ■■  r«Ur,Ur,  of  the  stream  surface,  is  composed  of  three 

parts:  is  the  inlet-outlet  of  the  flow;  is  the  interescting  line  of  the 

stream  surface  and  solid  wall  (blades,  rotor,  or  rotor  cover);  rp  is  the  artificial 
boundary  line.  should  satisfy 


yr 

-  g 


0lr#,  ■■  <Mrn  ■+*  (S, 


9m 


m. 

dm 


(2.20) 


In  the  equation,  (n^)  is  the  covariant  component  of  the  unit  vector  of  the  external 
normal  of  the  boundary. 

The  physical  components  W  ,  W  ,  and  W  (in  the  cylindrical  coordinate  system 

T  Ip  Z 

of  W)  can  be  derived  through  vP  ,  that  is 

By  using  equation  (2.18),  it  is  easy  to  derive  a  calculation  equation  (of  flow 
density)  bpW  _  *  (2.22) 

According  to  the  relationship  between  p  and  W,  there  is 
p  -  p.(c  -  (r  -  1 
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(2.23) 


In  the  equation,  c  —  l  —  (r  —  1)(«>V  —  2Mi)/2flJ;  Ps  and  ag  are,  respectively, 

the  stagnant  density  and  stagnant  sound  velocity  of  the  reference  point;  X  is  the 
degree  of  pre-rotation. 

Hence,  equations  (2.19)  and  (2.20)  are  nonlinear  boundary-value  problems. 

When  the  flow  is  subsonic,  then  these  are  nonlinear  elliptical  boundary-value 
problems. 


When  the  parameter  equations  are  determined  in  the  following  forms: 

r  —  r(x').  <P  —  <?(*')»  *  —  *(x“) 


then 


«<*  ”  rMr.,  +  rl<p^<p,t  +  x„x.# 
to1  —  rto(r.a«p.M‘  —  • 

If  the  stream  surface  is  a  rotating  surface,  let  z-z(m)  and  r=r(m)  be  the 

12 

parameter  equations  of  the  generatrix.  Then  by  assuming  x  =m  and  x  ^ , 

*i,  *  1»  Xu  “*  «n  ”  0,  *n  ™  “  r> 

<•  •  dr  /dm 

It  is  of  some  importance  that  the  blade  surface,  rotor  and  rotor  cover  in 
the  flow  channel  are  natural  stream  surfaces.  Hence,  these  stream  surface  equa¬ 
tions  are  known.  In  the  case  of  the  blade  surface,  if  its  parameter  equation  is 

1  2 

<j>=<t>(r,  z) .  At  that  time,  by  taking  x  =z  and  x  =r,  then 

«U  —  1  4-  rVi  ,  ■“  r JqMP.»  —  1  +  rV? 

m  —  1  +  r’Cvi  +  )>  «•’  “  —  rtotpj*/  • 

can  be  directly  solved  for  boundary-value  problems  (2.19)  and  (2.20)  in  order  to 
obtain  the  distribution  of  the  stream  field  at  the  surface  of  blade.  The  flow- 
field  distribution  at  the  blade  surface  is  the  main  basis  for  determining  the 
aerodynamic  characteristics  of  blade  types. 

III.  Approaching  a  Solution  of  Finite  Elements 

If  (*'*  x*)€0  and  V  is  a  Banach  space,  satisfying  H&0)(ZV<ZH'(Q), 
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Vmm  «|r,-0,  Jj  |rJ* 

then  corresponding  to  equations  (2.19)  and  (2.20),  the  simulated  dual-linear 
arbitrary  function  (over  VxV)  is 

B(^,  <£*)  -•  Wa^r-- (3.1) 
and  the  linear  arbitrary  function 

(Fg,  **>  -  jja  F<J>*Jx'Jx>  -  jr  ,*•*/  (3.2) 

We  assume  ifiH'iOi),  which  satisfies  ^1  —  <$|  ^  I  +  <»,  — I  — ~ 

lr,  I  >>i  I  Tf,  ft,  I  r„  Qn  J» 

Then,  the  weak  solution  of  the  boundary- value  problem  of  equations  (2.19)  and 
(2.20)  is:  we  solve  for  #  —  so  that 

B(*,4>m)-<Fgt*m>  V*«€F  (3.3) 

We  assume  that  V  £V  is  the  n-dimensional  subspace  and  its  fundamental  function 
system  is  The  Galerkin  approach  solution  of  equation  (3.3)  is  ex¬ 
pressed  as  .  *  . 

*■-  +  +  (3.4) 

Then  should  satisfy  the  algebraic  equation 

^v(«)**  —  fi(«)  "®  (33) 

In  the  equation 

A^e)  “  j  *M*,**k^*^*^  (3*) 

Fy(c)  —  —  \\a^^  ffoi+t**1*1*** +  |r# 


It  can  be  proved  that  the  solution  of  algebraic  equation  (3.5)  exists.  Hence, 
the  Galerkin  approach  solution  of  equation  (3.4)  also  exists.  It  can  be  proved 
that  the  solution  is  also  convergent. 

IV.  Calculation  Examples 

By  using  the  derived  differential  equations  (2.19)  and  (2.20)  for  boundary- 

value  problems,  solutions  are  obtained  by  using  the  finite-element  method  to 
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separately  calculate  flow  fields  of  a  T1102  blade  type  and  Gostelow  accurate* 
solution  blade  grid.  In  calculations,  the  finite-element  network  used  is  an 
arbitrary  quadrilateral  Con  a  curved  surface) ;  the  element  types  are  eight -nodal- 
point  quadrilateral  and  other  parameter  elements.  Compare  the  calculation  results 
with  that  of  the  T1102  blade  type  and  experimental  results  as  shown  in  Fig.  2; 
also  compare  the  results  with  that  of  the  Gostelow  blade  grid  and  theoretical 
solution  as  shown  in  Fig.  3.  Both  comparisons  match  to  each  other  relatively 
well. 


(C)  A  *mp  fit 

— ,  Wtllf  *■»  »*  9.9m  W*r- 
0.4H  91*  9.9H9  27. M 


Fig.  2.  Comparisons  of  experimental  values  with 
T1102  blade-type  plane  blade  grid  in  pressure 
distribution  along  the  blade  surface. 

Key:  (a)  Back  arc;  (b)  Internal  arc;  (c)  Calculated 
in  the  paper;  (d)  xx  experiments. 


Fig.  3.  Comparisons  with  accurate  solution 
of  velocity  distribution  at  blade  surface  of 

GOSTELOW  blade  grid:  _  accurate  solution; 

x  finite-element  solution  in  the  paper. 
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IRIDIUM  ELECTROSTATIC  PROBE  AND  ITS  APPLICATION  IN  A  HIGH  TEMPERATURE  WIND  TUNNEL 
Wang  Boyi 
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The  paper  describes  a  kind  of  electrostatic  probe  which 
can  be  used  to  measure  parameters  of  high  temperature  gas  flow 
for  such  parameters  as  electron  density  and  electron  temperature. 

The  paper  also  gives  theoretical  curves  of  probe  characteristics 
for  air  plasma.  This  iridium  electrostatic  probe  has  been 
operated  successfully  in  an  electric-arc-heated  wind  tunnel. 

I .  General  Descriptions 

In  experiments  and  studies  of  high  temperature  high  speed  aerodynamics,  as 
well  as  experiments  and  studies  of  an  aerodynamic  laser  and  electromagnetic 
fluid  dynamics,  the  electron  density  and  electron  temperature  of  a  high  tempera¬ 
ture  medium  are  very  important  parameters;  these  parameters  have  to  be  measured 
experimentally.  The  paper  introduces  a  type  of  electrostatic  probe  applied  in 
measuring  electron  density  and  electron  temperature  of  a  gas  stream  in  high 
temperature  wind  tunnels. 


'1 


The  electrostatic  probe  is  an  important  tool  for  plasma  diagnostics  [1] . 

The  measuring  and  testing  system  of  the  electrostatic  probe  is  simple,  with  a 
large  dynamic  range.  The  probe  can  measure  the  local  value  of  parameters.  The 
so-called  electrostatic  probe  is  actually  a  metallic  electrode.  The  probe  is 
placed  in  the  gas  to  be  measured;  the  external  power  source  provides  voltage 
to  obtain  a  certain  value  of  the  potential  relative  to  the  gas  to  be  measured. 

A  certain  number  of  electrons  and  ions  are  collected  from  the  gas  to  form  a 
current.  When  the  potential  of  the  probe  varies,  the  probe  current  will  vary 
according  to  a  fixed  rule.  After  testing,  the  ’’volt-ampere  characteristics"  of 
the  probe  are  obtained.  The  volt- ampere  characteristics  are  related  to  parameters 
such  as  density  and  temperature  of  charged  particles.  The  probe  theory  can  give 
the  quantitative  relationship  between  density  and  temperature  of  the  charged 
particles  in  the  gas. 


Fig.  1.  Electrostatic  probe  and  its  volt- 
ampere  characteristic  curve. 

Key:  (1)  Electrostatic  probe;  (2)  Plasma; 
(3)  Compensation  electrode. 
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In  the  severe  hot  environment  of  a  high  temperature  wind  tunnel,  first 
the  problem  of  survival  of  the  probe  in  the  high  temperature  stream  has  to  be 
solved  by  using  an  electrostatic  probe  for  measurement.  Hence,  we  developed 
an  iridium  electrostatic  probe,  and  it  has  been  successfully  operated  in  the 
H-4  electric-arc-heated  wind  tunnel  to  calibrate  electron  density  and  electron 
temperature  of  the  testing  stream  and  its  cross  section.  Basically,  the  calibra¬ 
tion  result  agrees  with  the  theoretical  calculation. 

II.  Structural  Design  of  Iridium  Electrostatic  Probe 

According  to  the  high- temperature  and  high-speed  characteristics  of  the 
electric-arc-heated  wind  tunnel,  the  electrostatic  probe  we  developed  is  composed 
of  three  parts:  central  electrode,  insulation  lining,  and  water-cooling  sleeve. 


Fig.  2.  Structural  diagram  of  iridium  electro¬ 
static  probe. 

Key:  (a)  Iridium  probing  electrode;  (b)  Boron 
nitride  insulation  lining;  (c)  Water-cooling 
sleeve;  (d)  Inlet  of  cooling  water;  (e)  Outlet 
of  cooling  water;  (f)  Insulation  sleeve;  (g) 
Lead-out  wire. 


The  central  electrode  includes  the  iridium  probe,  lead-out  wire,  heat- 
endurance  insulation  sleeve,  and  others.  Iridium  wire  has  the  advantages  of  high 
melting  temperature,  good  electrical  and  heat  conductivities,  high  energy  function, 
and  resistance  to  oxidation.  In  addition,  the  iridium  wire  is  hard  and  stiff, 
so  it  is  quite  adaptable  to  be  used  as  probes  for  operation  in  a  high-temperature 
air  environment.  As  proved  by  experiments,  this  iridium  electrostatic  probe 
can  be  successfully  and  continuously  operated  in  a  gas  stream  of  14.6  MJ/kg  of 
entropy  and  20  g/s  of  flow. 


Fig.  3.  Schematic  diagram  of  iridium 
electrostatic  probe  system  for  calibra¬ 
tion  of  electric-arc-heated  wind  tunnel. 
Key:  (a)  Direction  of  gas  stream;  (b) 
Nozzle;  (c)  Testing  sector;  (d)  Probe; 

(e)  Pressure  expander;  (f)  Control  system 
of  test  stand;  (g)  Water  cooling  system 
of  electrostatic  probe;  (h)  Power  feed 
system  of  electrostatic  probe;  (i)  Measure¬ 
ment  system  of  electrostatic  probe. 


A  high  pressure  (28  atmospheric  pressures)  water  pump  is  used  in  the  water¬ 
cooling  sleeve.  Besides  providing  protection  for  the  central  electrode,  the 
sleeve  is  also  used  as  a  compensation  electrode  of  the  probe. 


The  insulation  lining  can  ensure  electrical  insulation  between  the  external 
water  sleeve  and  the  iridium  electrode;  boron  nitride  is  used  for  the  lining,  as 


the  material  is  high  in  insulation,  bum  etching  and  heat  conduction  character¬ 
istics.  Therefore,  some  cooling  effect  can  be  provided  for  the  iridium  wire  by 
these  characteristics.  A  high-temperature  seal  is  used  between  the  central 
electrode  and  boron  nitride. 

The  entire  probe  is  placed  on  a  test  stand  of  a  three-dimensional  wind 
tunnel . 

III.  Theory  and  Data  Analysis  of  Electrostatic  Probe 


The  theoretical  analysis  of  the  electrostatic  probe  is  relatively  complex. 
In  particular,  in  calibration  of  the  gas  stream  in  a  wind  tunnel,  the  flow 
effect  and  electrostatic  effect  of  the  plasma  should  be  coupled;  thus,  the 
theoretical  description  is  even  more  complex.  In  order  to  avoid  considering 
the  flow  effect,  we  used  a  cylindrical  probe  with  its  axial  line  parallel  to  the 
direction  of  the  gas  stream  velocity.  In  addition,  the  probe  diameter  should 
be  sufficiently  small  in  order  to  satisfy  the  theoretical  requirements  of  the 
electrostatic  probe  of  free  molecules: 

DtR,<l 

In  the  expression,  D  is  the  length  of  the  Durban  screen;  is  the  probe  radius; 
and  X  is  the  average  free  range  of  charged  particles  in  the  gas  stream. 


High  temperature  gas  is  a  partially  ionized  plasma;  in  the  plasma  the  densi¬ 
ties  of  electrons  and  ions  are  equal  to  each  other,  n  ~n.  (n  and  n.  are, 

6  1  6  X 

respectively,  densities  of  electrons  and  ions) .  Hence,  the  probe  current  is 


composed  of  electron  current  and  ion  current:  Ip=Ie-I^. 


In  the  equation. 


I 

P 


is 


the  probe  current;  I  is  electron  current;  and  I,  is  ion  current.  We  can  see 

©  A 


from  Fig.  1  that  when  the  probe  potential  is  lower  than  the  space  potential 
VQ  of  the  plasma  (generally,  it  is  called  the  negative  probe) ,  the  probe  will 


attract  ions  but  reject  electrons.  The  higher  the  negative  value  of  the  potential. 


the  more  are  electrons  rejected.  Hence,  the  probe  current  drops  rapidly.  At  the 


floating  potential  V^,  the  electron  current  and  ion  current  are  equal  in  value 
and  the  probe  current  is  zero.  When  the  probe  potential  is  considerably  lower 
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than  the  floating  potential,  the  probe  collects  a Isos t  all  the  ion  current. 
Besides,  the  current  varies  very  slowly  with  the  potential;  generally,  such  a 
current  is  called  the  saturated  ion  current. 


An  accurate  probe  theory  [2]  gives  the  expression  for  the  ion  current: 

(££•)*  /♦<*»«.«)  (I) 

In  the  equation,  is  the  density  of  the  ion  current  of  the  probe;  e  is  the 
electron  charge;  T  is  the  electron  temperature;  m.  is  the  ion  quality;  k  is 
Boltzmann's  constant;  I+  is  the  current  density  of  dimensionless  ions;  and  Xpt 
£  and  e  are  dimensionless  variables.  Their  definitions  are: 


In  the  equations,  is  the  ion  temperature.  In  paper  [2],  a  numerical  method 
is  used  to  obtain  the  dimensionless  function  I+  as  the  calculation  results  of 
the  functions  xp»  €  and  e.  Starting  from  this  result,  we  derived  theoretical 
curves  (refer  to  Fig.  4)  adaptable  to  air  plasma;  the  abscissa  is  f*f+(jr#  — 10), 

and  the  ordinate  is  /•*.(*,  —  10).  ’  Here,  which  is  the  dimension- 

less  floating  potential.  I+(x£-10)  indicates  that  the  potential  is  lower  than  the 
floating  potential  by  a  dimensionless  value  of  densities  of  ions  and  current  at 
10(kT  /e)  volts.  When  the  electron  temperature  T  is  known,  in  the  volt-ampere 
characteristics  of  the  probe  measured  by  experiments,  it  is  not  difficult  to 
determine  the  position  of  (X£-10);  the  value  of  I ^ (x £-10)  can  be  found  accordingly. 

Thus,  it  is  required  to  first  determine  electron  temperature  from  plotted 
curves  as  shown  by  the  two  following  methods: 


1.  The  inclination  of  a  semi- logarithmic  curve  of  electron  current  of  the  negative 


probe: 


4L- 


2.  The  inclination  of  the  probe  characteristic  curve  at  the  floating  potential 


il>. 


[iit  - 1  i£<|] 

ldV,  dV,J 


\wr-rt 


(3) 


Electron  temperature  values  can  be  simply  determined  by  using  a  graphic  method 
(Fig.  S). 


Fig.  4.  Electrostatic  probe: 
theoretical  curves  in  air  plasma. 


Fig.  5.  Graphic  method  of  determining 
electron  temperature. 

Key:  (a)  Inclination  of  volt-ampere 
characteristic  of  probe  at  floating 
potential;  (b)  Extension  line  of 
saturated  ion  current  of  the  probe; 

(c)  Electron  temperature  value  determined 
by  graphic  method. 


Given  a  known  electron  temperature,  we  can  utilize  the  measured  probe 
characteristics  and  the  given  theoretical  curves  to  determine  the  density  of 
electrons  or  ions.  We  use  the  following  steps: 


1.  From  the  measured  ion-current  curve  of  the  probe,  I^(Xf-lO)  can  be  determined. 

2.  We  calculate  the  density  of  the  ion  current  at  potential  (x£-10) 

/,<*/  -  io)  -  (4) 

In  the  equation,  Ap-2irRp£.  is  the  surface  area  of  the  probe;  i  is  the  length  of 
the  probe. 

2 

3.  We  calculate  the  corresponding  value  of  5  I+(x£-10) 


?!+(X,  -  10)  -  1.563  X  10**1 

■  r 


In  the  equation,  the  unit  of  Rp  is  M;  the  unit  of  is  A/M  ;  and  the  unit  of 

kT  /e  is  V. 
e 

4.  From  the  theoretical  curve  (Fig.  4)  of  the  probe,  the  I+fyf-10)  value  of 
52I+Ocf-10)  is  determined. 

5.  We  calculate  the  corresponding  value  of  the  ion  density 

_  _  1, Or#  - 10) 


-*** 


In  the  equation,  the  unit  of  ^  is  [1/M3];  other  units  are  the  same  as  above. 

For  plasma,  n  *n.;  therefore,  the  electron  density  is  also  determined. 

0  i 
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IV.  Calibration  Result  of  Electron  Density  in  Gas  Stream  During  Wind  Tunnel 
Testing  Sector 


By  using  the  iridium  electrostatic  probe,  calibration  tests  were  conducted 
in  the  H-4  electric-arc-heated  wind  tunnel.  The  basic  characteristic  parameters 
of  the  H-4  wind  tunnel  are  as  follows: 


Electric  power  source:  1500  kW;  Work  medium:  air;  Gas  flow:  15-40  g/s;  Gas 

5  6  2 

pressure  in  electric-arc-heated  chamber:  3x10  --10  N/m  ,  and  total  entropy  of 
gas:  8-18  MJ/kg.  In  addition,  the  throat  diameter  of  the  nozzle  is  8.1  mm;  the 

exit  diameter  is  220  mm;  the  total  pressure  in  the  testing  sector  during  measure- 

2  2  2 
ment  is  approximately  3x10  N/m  ;  and  the  static  pressure  is  approximately  3  N/m  . 


According  to  the  operation  parameters  of  the  wind  tunnel  mentioned  above, 

we  calculated  the  theoretical  values  of  the  electron  density  and  electron 

temperature  of  the  test  gas  stream  (during  estimation  and  calculation,  the 

5  2 

pressure  of  the  electric-arc-heated  wind  tunnel  is  assumed  to  be  6.87x10  N/m  , 
that  is  6.8  atmospheric  pressures).  Results  are  shown  in  Figs.  6  and  7. 


Fig.  6.  Estimation-calculation  of  Fig.  7.  Estimation-calculation  of  electron 
electron  density  of  test  gas  stream  temperature  range  of  test  gas  stream  of  H-4 
in  H-4  electric-arc-heated  wind  tunnel,  heated  chamber. 

Key:  (a)  Electron  density;  (b)  Entropy  Key:  (a)  Electron  temperature;  (b)  Entropy 
value  of  wind  tunnel;  (c)  Frozen  flow;  value  of  wind  tunnel;  (c)  Total  atmospheric 
(d)  Balanced  flow;  (e)  Total  atmo-  pressure  in  electric-arc-heated  chamber, 
spheric  pressure  in  the  heated  chamber.  50 


Calibration  testing  proceeds  in  two  stages.  The  first  stage  emphasizes  the 
measurement  of  local  values  of  parameters  such  as  electron  density.  At  three 
different  cross  sections  of  the  gas  stream,  testing  points  in  all  were  tested. 

The  entire  volt-ampere  characteristic  is  obtained  by  scanning  with  a  voltage 
probe.  After  compilation,  typical  measurement- testing  results  are  plotted  in 
Figs.  8  and  9.  We  know  from  these  measurement  results  that  the  test  gas  stream 
of  the  H-4  electric-axc-heated  wind  tunnel  approaches  a  frozen  stream.  The 
measured  electron  density  basically  matches  the  degree  of  magnitude  of  the 
theoretical  estimate  calculation.  However,  the  measured  value  of  the  electron 
temperature  is  higher  than  the  estimate  calculation.  This  is  possibly  because 
the  graphic  solution  method  is  not  too  accurate.  At  the  same  time,  as  revealed 
by  the  measurement  result  of  the  electron  temperature,  basically  the  electron 
temperature  of  the  test  gas  current  is  homogeneous  because  the  measured  values 
at  various  measurement  points  are  randomly  distributed  within  the  error  range  of 
measurement.  This  result  matches  the  physical  characteristics.  Since  the  electron 
mass  is  considerably  smaller  than  the  mass  of  ions  and  neutral  molecules  in  the 
gas  stream,  during  their  collision  process  there  is  little  exchange  of  eneTgy 
between  the  electrons  and  these  particles.  Therefore,  the  electrons  are  unable 
to  maintain  energy  equilibrium  with  other  heavy  particles  and  still  maintain  the 
energy  state  in  front  of  the  throttle. 
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Fig.  8.  Calibration  result  of  electron  density  of 
test  gas  stream  in  H-4  electric-arc-heated  wind  tunnel. 
(The  cross  section  of  the  gas  stream  is  8  cm  from  the 
exit  of  the  nozzle;  the  entropy  value  of  the  wind  tun¬ 
nel  is  2870-3250  calories  per  gram.) 

(Key  on  the  following  51  page] 


[Key  of  Fig.  8  on  the  preceding  page] 

Key:  (a)  Center  line  of  testing  sector  of  wind  tunnel; 

(b)  Distances  (cm)  of  various  measurement  points  from 

the  center  line;  (c)  At  edge  of  nozzle;  (d)  Electron 

density  n  (cm_3) . 
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Fig.  9.  Calibration  and  measurement  results  of  electron 
temperature  of  test  gas  stream  in  H-4  electric-arc-heated 
wind  tunnel:  The  results  were  totalled  from  three  cross 
sections  0,  4  and  8  cm  from  the  exit  of  the  nozzle;  the 
entropy  value  of  the  wind  tunnel  is  2400-3900  calories  per 
gram. 

Key:  (a)  Center  line  of  testing  sector  of  wind  tunnel; 

(b)  Distances  (cm)  of  various  measurements  from  center  line; 

(c)  At  edge  of  nozzle;  (d)  Electron  temperature  T  (°K) . 
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The  calibration-measurement  testing  in  the  second  stage  emphasized  the 
distribution  of  measuring  the  electron-density  cross  section.  Since  it  can  be 
considered  that  within  the  same  cross  section  the  electron  temperature  is 
homogeneous,  the  variation  in  ion  current  represents  the  variation  in  electron 
density.  Thus,  the  probe  deviates  in  fixation  (voltage  fixed  at  -12  V)  and 
the  three-dimensional  support  can  continuously  move  across  the  gas  stream.  In 
such  a  way,  the  value  distribution  of  the  saturated  ion  current  measured  is 
the  distribution  shape  of  the  cross  section  of  the  electron  density.  Figure  10 
shows  a  typical  measurement  curve. 
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Fig.  10.  Measurement  curve  of  distribution  of  saturated 
ion-current  cross  section  in  test  gas  stream  of  H-4 
electric-arc-heated  wind  tunnel:  the  curve  shows  the 
cross  sectional  result  15  cm  from  the  nozzle  exit. 

Key:  (a)  Center  line  of  testing  sector  of  wind  tunnel; 

(b)  Diametral  distance  (cm) ;  (c)  At  edge  of  nozzle; 

(d)  Saturated  ion  current  1^  (mA)  of  probe. 

V.  Conclusion 

The  testing  of  calibration  and  measurement  indicates  that  the  iridium 
electrostatic  probe  can  be  operated  in  a  high-temperature  wind  tunnel;  the 
characteristics  of  the  probe  meet  the  testing  requirements. 

The  electrostatic  probe  can  be  used  to  measure  the  electron  density,  electron 
temperature,  and  the  distribution  of  the  cross  section  of  a  gas  stream  tested  in 
a  wind  tunnel.  At  the  same  time,  results  of  probe  measurement  also  provide  a 
basis  for  evaluating  the  degree  of  nonequilibrium  of  a  gas  stream  in  a  wind 
tunnel  because  the  difference  in  degree  of  nonequilibrium  of  the  gas  stream  in 
the  wind  tunnel  can  apparently  be  revealed  in  the  magnitude  of  the  electron 
density  in  the  gas  stream.  By  assuming  the  degree  of  nonequilibrium  of  the  gas 
stream  in  the  wind  tunnel  and  calculation  of  electron  density  based  on  theoretical 
calculation,  and  comparing  with  the  measured  electron  density,  the  degree  of  non¬ 
equilibrium  of  the  gas  stream  can  be  determined. 
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